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Abstract. Explosive astrophysical systems, such as supernovae or compact star 
binary mergers, provide conditions where strange quark matter can appear. The high 
degree of isospin asymmetry and temperatures of several MeV in such systems may 
cause a transition to the quark phase already around saturation density. Observable 
signals from the appearance of quark matter can be predicted and studied in 
astrophysical simulations. As input in such simulations, an equation of state with 
an integrated quark matter phase transition for a large temperature, density and 
proton fraction range is required. Additionally, restrictions from heavy ion data and 
pulsar observation must be considered. In this work we present such an approach. We 
implement a quark matter phase transition in a hadronic equation of state widely used 
for astrophysical simulations and discuss its compatibility with heavy ion collisions and 
pulsar data. Furthermore, we review the recently studied implications of the QCD 
phase transition during the early post-bounce evolution of core-collapse supernovae 
and introduce the effects from strong interactions to increase the maximum mass of 
hybrid stars. In the MIT bag model, together with the strange quark mass and the bag 
constant, the strong coupling constant a s provides a parameter to set the beginning 
and extension of the quark phase and with this the mass and radius of hybrid stars. 



Strange quark matter in explosive astrophysical systems 



2 



Figure 1. The phase di- 
agram for the quark mat- 
ter phase transition for 
SN matter (solid lines) 
and HI collisions (dashed 
lines). Thin lines de- 
note the beginning of the 
mixed phase, thick lines 
mark the onset of the pure 
quark phase. 

1. Introduction: 

The future FAIR facility at GSI, Darmstadt, will explore the equation of state (EoS) of 
strongly interacting matter for intermediate temperatures T and high baryon densities 
nb around isospin symmetry, that is for proton fractions Y p ~ 0.5. Supernovae (SNe) 
and binary mergers hold environments with similar conditions for T and but with 
Y p < 0.3. As will be discussed in the scope of this article, core-collapse SNe with matter 
at a low value of Y p and dynamical timescales in the range of ms, provide conditions 
suitable for a phase transition to strange quark matter. Such a scenario was recently 
studied in p] applying the MIT bag approach for the EoS of quark matter and using 
low critical densities for its onset. Simulations with different progenitor models and two 
different bag constants led to SN explosions accompanied by a significant neutrino burst 
which can be observed by present and future neutrino detectors [2]. In the following, 
we will introduce in more detail the hybrid EoS used in the above work and analyze 
its influence on the dynamics of the PNS evolution. We will discuss the compatibility 
with heavy ion (HI) data and pulsar observations. Furthermore we will include first 
order corrections from the strong interaction constant a s and study its influence on the 
maximum mass of the cold hybrid star configurations. 

2. Initial setup for the equation of state 

For the hadronic part of the quark-hadron EoS we use the relativistic mean field 
approach by (3], while quark matter is described by the MIT bag model. Due to their 
small current mass, the up and down quarks are treated as massless, while for the strange 
quark we chose m s = lOOMeV which is well within the limits set by the Particle Data 
Group jl]. If quark masses are fixed and no corrections from the strong coupling constant 
a s are included, the critical densities for the phase transition are directly given by the bag 
parameter B for which we applied two values, B 1 ^ = 165MeV and B 1 ^ = 162MeV. For 
the construction of the phase transition we choose the Gibbs approach, where a mixed 
phase of quarks and hadrons is present (5J E] . Figure [l] shows two phase diagrams with 
the pure hadronic and quark phases and the mixed phase for B 1 ^ = 165MeV. As can be 
seen, the onset of strange quark matter (denoted as uds) for a proton fraction of Y p = 0.3 
happens already around saturation density no- For B 1 ^ 4 = 162MeV the critical densities 
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are even smaller. However, these values do not contradict with results from HI collisions 
due to two main reasons. First, supernova dynamics happen on timescales of ms, whereas 
weak processes operate within 1CT 6 — 10~ 8 s and have therefore enough time to produce 
strangeness. Consequently, phase transitions can be considered from hadronic to three 
flavor quark matter in weak equilibrium [7] . For HI collisions, dynamical timescales are 
much shorter, of the order of 10~ 23 s and therefore not long enough for strangeness to 
be produced and equilibrated by weak interactions. Consequently, for such systems, it 
seems to be more appropriate to consider a phase transition from hadronic to quark 
matter composed only of up and down quarks. The higher the number of quark flavours 
and therefore the number of degrees of freedom, the softer is the EoS in the mixed phase 
and the lower is the critical density for its onset. The second main difference for SN 
and HI environments is the proton fraction in the two systems, being Y p < 0.3 for the 
first and Y p ~ 0.5 in the second case. Due to the symmetry energy of hadronic matter, 
its isospin symmetric state is energetically favored. For a proton fraction Y p < 0.5, 
the energy of hadronic matter is higher and the additional asymmetry pressure stiffens 
the EoS. This stiffness results in an earlier onset of the mixed phase with its softer 
EoS. Consequently, a low value of Y p leads to smaller critical densities than for isospin 
symmetric matter. Figure [T] shows two phase diagrams, for SN environments and HI 
collisions, for B l l 4 = 165MeV, illustrating that a low onset of quark matter in SN 
environments is compatible with a high critical density in HI collisions, which, for the 
chosen B and small T, is up to 5n . However, the exact location of the critical density 
for different T varies in dependence of the models for the quark and hadron EoSs, or 
the inclusion of finite size effects. 



The softening in the mixed phase, as seen in figure 2(a) , is caused by the growing 



quark fraction x- Figure 2(b) shows the fractions of positive charge Yo in quark and 
hadronic matter in the pure and mixed phases. In the quark phase Yq is given by 
(2/3 n u — 1/3 nd — 1/3 n s ) /rib, whereas n u , and n s are the up, down and strange 
quark number densities. For the hadronic phase, the charge fraction corresponds to the 
proton fraction Y p . As shown in figure 2(b) , Yq in the quark phase can be very low 
and even negative. Therefore, with increasing x, the charge fraction in the quark phase 
can compensate Yc = Y p of hadronic matter, and the latter can consequently approach 
isospin symmetry towards the end of the mixed phase. At this point, due to the soft EoS 
in the isospin symmetric hadronic phase and the large number of degrees of freedom, 
the mixed phase EoS is very soft. However, the vanishing of hadronic degrees of freedom 
causes a significant stiffening when the pure quark phase sets in. Nevertheless, in the 
simple MIT bag model, the EoS of the pure quark phase is still much softer than the 
one for hadronic matter giving low maximum masses for hybrid stars. 
Up to now, the highest precisely measured mass is the one for the Hulse- Taylor pulsar 
with 1.4414 ± 0.0002 solar masses M Q . A new candidate might be the recently studied 
J1903+0327, a millisecond pulsar with a main sequence star companion. Due to the 
large eccentricity of the binary system, the advance of periastron can be measured giving 
a value of 1.67±O.O1M for the mass of J1903+0327 |8J. The mass-radius relations in 
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Figure 2. (a) Hybrid EoS for B 1 / 4 = 165MeV and a s = 0, together with the EoSs for 
the pure hadronic and quark phases at Y p = 0.3 and T = 15MeV; (b) Charge fractions 
in quark and hadronic matter in the mixed phase for a total Yq = 0.3 and temperature 
T = 15MeV using B 1 / 4 = 165MeV and a s = 0. 




Figure 3. (a) The inclusion of first order corrections from the strong interaction 
constant a s for quarks in the MIT bag model leads to an increase in the critical 
density. This can be compensated by reducing the value of B. The reduction in B 
and inclusion of a s results in a higher pressure in the mixed and quark phases and 
therefore a higher hybrid star maximum mass (b). 



figure 3(b) show that while hybrid stars for B 1 ^ = 165MeV and B l l A = 162MeV are 
above the Hulse- Taylor constraint of 1.44M , their maximum masses are smaller than 
1.67±0.01M Q . A possibility to increase the maximum masses of hybrid stars within the 
MIT bag model is the inclusion of first order corrections from the strong interaction 
constant a s [9]. Figure 3(a) shows on the example of B 1 ^ 



165MeV how the inclusion 
of a s corrections shifts the critical density to higher values, increasing the pressure in 
the mixed and the pure quark phases. A low density for the onset of quark matter can 
be obtained again by decreasing the value for B. However, the new parameter set for 
B and a s leads to a stiffer EoS in the mixed and quark phases which results in higher 
maximum masses of hybrid stars. This is shown in figure 3(b) where the parameter set 
of B 1 / 4 = 155MeV with a s = 0.3 leads to a hybrid star maximum mass of ~ 1.67M 
and, at the same time, has a similar critical density for the onset of quark matter as 
B 1 / 4 = 165MeV with a, = 0. 
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3. Quark matter in supernovae 

We performed core-collapse simulations of low and intermediate mass Fe-core progenitors 
in spherical symmetry with two different bag constants B 1 ^ = 165MeV and B 1 / 4 = 
162MeV for the quark matter EoS. Our numerical model is based on general relativistic 
radiation hydrodynamics and three flavor Boltzmann neutrino transport (for details see 
[TO] and references therein). The conditions for the appearance of quark matter, i.e. 
the beginning of the mixed phase, are already obtained at the Fe-core bounce at central 
densities close to and slightly above n . However, the produced small quark fraction 
initially does not influence the SN dynamics and the evolution proceeds like in a normal 
core collapse supernova. A hydrodynamic shock wave forms, travels outwards and looses 
energy due to the disintegration of infalling heavy nuclei and production of neutrinos. 
The latter become observable in the neutrino spectra as a neutrino burst dominated 
by electron neutrinos as the shock wave propagates across the neutrinospheres (i.e. 
the neutrino energy and flavor dependent spheres of last scattering). These energy 
losses turn the expanding shock quickly into a standing accretion shock already ~ 5ms 
after bounce. As discussed previously, the reason for the unchanged dynamics are the 
relative stiffness and similarity of the mixed phase EoS to the hadronic one at small 
quark fractions and Y p close to isospin symmetry. However, as matter continues to 
be accreted on the surface of the PNS, the density and temperature in its interior 
increases and a growing volume enters the mixed phase. The quark matter fraction in 
the mixed phase rises and the PNS interior moves up to softer regions of the EoS, where 
it becomes gravitational unstable. A contraction proceeds into a collapse till pure quark 
matter is reached and the EoS stiffens again due to the disappearance of the additional 
hadronic degrees of freedom. The collapse halts and a second shock front forms at 
the phase boundary between the mixed and hadronic phases. This second shock front 
moves outwards and turns into a shock wave when it reaches the PNS surface where 
the density drops over several orders of magnitude. Hereby, the decrease in density 
accelerates the shock wave to velocities of the order of the vacuum speed of light. Shock 
heating of infalling hadronic matter leads to a lift of degeneracy and an increase in its 
proton fraction accompanied by the production of anti-neutrinos. As soon as the shock 
wave propagates over the neutrinospheres, a second neutrino burst, dominated by anti- 
neutrinos is released. The delay of this second burst after the first deleptonization burst 
contains correlated information about the progenitor model, the hadronic and quark 
EoSs, and the quark-hadron phase transition. For more details, see [TJ. 

4. Conclusions 

Because of the different proton fractions of matter in terrestrial and astrophysical 
laboratories, such as the future FAIR facility at GSI and supernovae or compact star 
mergers, on the one hand and their similarities in T and n& on the other, the study of 
heavy ion collisions and explosive astrophysical scenarios can complement each other in 
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probing the phase diagram of strongly interacting matter, also in regard to the phase 
transition from hadronic to quark matter. However, the study of possible observable 
signals and impacts of quark matter in astrophysical systems requires hydrodynamical 
simulations with an input of an appropriate quark-hadron equation of state. In this 
article we present such an approach where a quark matter phase transition has been 
implemented in a hadronic equation of state for a large range of temperatures, proton 
fractions and densities. Applying the latter to simulations of core-collapse supernovae, 
we find that a quark matter phase transition can cause the formation of a second shock 
wave which leads to the explosion of the star, accompanied by a second neutrino burst. 
The latter is dominated by anti-neutrinos, which can be observed by future and present 
neutrino detectors. If found, the second neutrino peak can give correlated information 
about the progenitor mass and the critical density for the onset of quark matter. 

Acknowledgement The project was funded by the Swiss National Science Foundation 
grant, no. PP00P2- 124879/1 and 200020-122287 and the Helmholtz Research School 
for Quark Matter Studies, the Italian National Institute for Nuclear Physics, the 
Graduate Program for Hadron and Ion Research (PG-HIR), the Alliance Program of the 
Helmholtz Association (HA216/EMMI) and the DFG through the Heidelberg Graduate 
School of Fundamental Physics. The work of G. P. is supported by the Deutsche 
Forschungsgemeinschaft (DFG) under Grant No. PA 1780/2-1. The authors are 
additionally supported by CompStar, a research networking program of the European 
Science Foundation, and the Scopes project funded by the Swiss National Science 
Foundation grant, no. IB7320-110996/1. 

[1] I. Sagcrt, T. Fischer, M. Hempel, G. Pagliara, J. Schaffner-Biclich, A. Mezzacappa, F.-K. 
Thielemann, and M. Liebendorfer. Phys. Rev. Lett., 102, February 2009. 

[2] B. Dasgupta, T. Fischer, S. Horiuchi, M. Liebendoerfer, A. Mirizzi, I. Sagert, and J. Schaffner- 
Biclich. Detecting the QCD phase transition in the next Galactic supernova neutrino burst. 
ArXiv e-prints, December 2009. 

[3] H. Shcn, H. Toki, K. Oyamatsu, and K. Sumiyoshi. Nucl. Phys., A637:435-450, 1998. 

[4] C. Amsler. Particle Data Group. Phys. Lett. B, 667, 2009. 

[5] M. Hempel, G. Pagliara, and J. Schaffner-Bielich. Phys. Rev. D, 80, December 2009. 
[6] N. K. Glcndenning. Phys. Rev. D, 46, August 1992. 

[7] B. W. Mintz, E. S. Fraga, G. Pagliara, and J. Schaffner-Bielich. Nucleation of quark matter in 

protoneutron star matter. ArXiv e-prints, October 2009. 
[8] P. C. C. Freire. Eccentric Binary Millisecond Pulsars. ArXiv e-prints, July 2009. 
[9] E. Farhi and R. L. Jaffc. Phys. Rev. D, 30, December 1984. 
[10] M. Liebendorfer, O. E. B. Messer, A. Mezzacappa, S. W. Bruenn, C. Y. Cardall, and F.-K. 

Thielemann. Astrophys. J., SuppL, 150, January 2004. 



